Phase Diagram of spin 1/2 XXZ Model With Dzyaloshinskii-Moriya Interaction 



OO 
O 

o 

(N 

o 

Q 

o 



CO 

-4-3 
o3 



i 

C 

O 
o 



> 

(N 

\o 

OO 

(N 
OO 

o 



R. Jafari 1 ' 2 and A. Langari 3 

'institute for Advanced Studies in Basic Sciences, Zanjan 45195-1159, Iran 
^School of physics, IPM (Institute for Studies in Theoretical Physics and Mathematics), 
P. 0. Box: 19395-5531, Tehran, Iran 
s Physics Department, Sharif University of Technology, Tehran 11155-9161, TrarQ 

(Dated: December 10, 2008) 

We have studied the phase diagram of the one dimensional XXZ model with Dzyaloshinskii-Moriya 
(DM) interaction. We have applied the quantum renormalization group (QRG) approach to get the 
stable fixed points, critical point and the scaling of coupling constants. This model has three phases, 
ferromagnetic, spin-fluid and Neel phases which are separated by a critical line which depends on 
the DM coupling constant. We have shown that the staggered magnetization is the order parameter 
of the system and investigated the influence of DM interaction on the chiral ordering as a helical 
magnetic order. 
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INTRODUCTION 



Quantum phase transition has been one of the most 
interesting topic in the area of strongly correlated sys- 
tems in the last decade. It is a phase transition at 
zero temperature where the quantum fluctuations play 
the dominant role 1 . Suppression of the thermal fluctu- 
ations at zero temperature introduces the ground state 
as the representative of the system. The properties of 
the ground state may be changed drastically shown as a 
non-analytic behavior of a physical quantity by reaching 
the quantum critical point (QCP). This can be done by 
tuning a parameter in the Hamiltonian, for instance the 
magnetic field or the amount of disorder. The ground 
state of a typical quantum many body systems consist 
of a superposition of a huge number of product states. 
Understanding this structure is equivalent to establish- 
ing how subsystems are interrelated, which in turn is 
what determines many of the relevant properties of the 
system. In Mott insulators the Heisenberg interaction 
is in most cases the dominant source of coupling be- 
tween local moments, and most theoretical investiga- 
tions are based on modeling in which only this type of 
interaction is included. Recently some novel magnetic 
properties in antiferromagnetic (AF) systems were dis- 
covered in the variety of quasi-one dimensional materi- 
als that are known to belong to an antisymmetric in- 
teraction of the form D .(Si x Sj) which is known as 
the Dzyaloshinskii-Moriya (DM)interaction. The rel- 
evance of antisymmetric superexchange interactions in 
spin Hamiltonian which leads to either a week ferro- 
magnetic (F) or helical magnetic distortion in quantum 
AF systems, has been introduced phenomenologically by 
Dzyaloshinskii 2 . A microscopic model of antisymmet- 
ric exchange interaction was first proposed by Moriya 3 
which showed that such interactions arise naturally in 
perturbation theory due to the spin-orbit coupling in 
magnetic systems with low symmetry and is essentially 
an extension of the Anderson superexchange mechanism 4 



that shows for spin-flip hopping of electrons. Since it 
(DM interaction) breaks the fundamental SU(2) sym- 
metry of the Heisenberg interactions, it is at the ori- 
gin of many deviations from pure heisenberg behavior 
such as canting- or small gaps&i&^iiii. A number of 
AF systems expected to be described by DM interac- 
tion, such as Cu(C 6 D 5 COO) 2 3D 2 0&d±, YbiAs^^., 
BaCu 2 Si 2 Or^; a - Fe 2 3 , LaMnO^ and K 2 V 3 Og^., 
which exhibit unusual and interesting magnetic proper- 
ties in the presence of quantum fluctuations and/or ap- 
plied magnetic fiel d 16 : 18 ' 19 . Also belonging to the class 
of DM antifcrromagnets is La 2 CuO±, which is a parent 
compound of high-temperature superconductors 2 - . This 
has stimulated extensive investigation on the physical 
properties of the DM interaction. However, This interac- 
tion is rather difficult to Handel analytically, which has 
brought much uncertainty in the interpretation of exper- 
imental data and has limited our understanding of many 
interesting quantum phenomena of low-dimensional mag- 
netic materials. 



In the present paper, we have considered the one 
dimensional XXZ model with DM interaction by im- 
plementing the quantum renormalization group (QRG) 
method. In the next section the QRG approach will be 
explained and the renormalization of coupling constant 
are obtained. In section (III), we will obtain the phase 
diagram, fixed points, critical points and calculate the 
staggered magnetization as the order parameter of the 
underling quantum phase transition. We will also intro- 
duce the chiral order as an ordering which is produced 
by DM interaction. The exponent which shows the di- 
vergence of correlation function close to the critical point 
(v), the dynamical exponent (z) and the exponent which 
shows the vanishing of staggered magnetization near the 
critical point (/?) will also be calculated. In Sec (IV), 
discussion concludes the paper. 
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FIG. 1: (color online)The decomposition of chain into three 
site blocks Hamiltonian (H B ) and inter- block Hamiltonian 
(# ss ). 



II. QUANTUM RENORMALIZATION GROUP 

The main idea of the RG method is the elimination 
or the thinning of the degrees of freedom caried out step 
by step in an iteration procedure. Here, we used the well 
known Kadanoff block method as it is both well suited to 
perform analytical calculations in the lattice models and 
is conceptually straight-forward to be extended to the 
higher dimension s 21 ' 22 ! 23 ' 24 . In the Kadanoff's method, 
the lattice is divided into blocks in which the Hamilto- 
nian can be exactly diagonalized. Selecting a number of 
low-lying eigenstates of the blocks the full Hamiltonian 
is projected onto these eigenstates and an the effective 
(renormalized) Hamiltonian is obtained. 

The Hamiltonian of XXZ model with DM interaction 
in the z direction on a periodic chain of N sites is 



j N 

H(J, A) = - £ [a?af +1 + o*o* +1 + Aa*af +1 (1) 

i 

where the J is exchange constant, D is the strength 
of z component of DM interaction and the easy-axis 
anisotropy defined by A which can be positive and neg- 
ative. The positive and negative J corresponds to the 
antiferromagnetic and ferromagnetic (F) cases, respec- 
tively, uf refers to the cv-component of the Pauli matrix 
at site i. By implement n rotation around z axis on 
odd or even sites, the AF case of Hamiltonian (J > 0) 
is mapped on the F case (J < 0) with opposite sign of 
anisotropy, 



H(J, A) = - £ + afaf +1 - Aof a* +1 

i 



(2) 



J > 0. 



So we can restrict ourselves to AF case ( J > 0) with 
D > and arbitrary anisotropy (A < and A > 0) 
without loss of generality. 

The effective Hamiltonian to the first order RG ap- 
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FIG. 2: (color online) Phase diagram of the XXZ model with 
DM interaction. The long dashed line is the critical line which 
separates SF-Neel phases and is characterized by A c = (1 + 
D 2 ) 1//2 . The dashed dot line shows the level crossing which 
separate SF-F phases is given by A = — (1 + D 2 ) 1 ^ 2 . Arrows 
show the running of coupling constant under RG iteration. 
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Hl ff = P H BB P . 



We consider a three-site block procedure defined in 
Fig.©. The block Hamiltonian (H B = J2 h f)> its 
eigenstates and eigenvalues are given in Appendix A. 
The three-site block Hamiltonian has four doubly de- 
generate eigenvalues (see Appendix A). Pq is the pro- 
jection operator of the ground state subspace defined by 
{Pa = I 1t)(H + \ Wo I), Where |V>o> and |^) are the 
doubly degenerate ground states, | ft) and | JJ.) are the re- 
named base kets in the effective Hilbert space. For each 
block we keep two states (|V>o) an d IV'o)) to define the ef- 
fective (new) site. Thus, the effective site can be consid- 
ered as having a spin 1 /2. Due to the level crossing which 
occurs for the eigenstates of the block Hamiltonian, the 
projection operator (Pq) can be different depending on 
the coupling constants. Therefore, we must specify dif- 
ferent regions with the corresponding ground states. As 
Fig. ([5]) shows there are two regions with different eigen- 
states which are separated by A < — \/l + D 2 where a 
level crossing occurs. In region (A) the ground state is 
the doubly-degenerate ferromagnetic state ^3) and ^3) 
while in region (B) |?/>o) and \ip'o) are the degenerate 
ground states. At the level crossing (A = —yl + D 2 ) 
the ground state is 4- fold degenerate ( | ^3 ) , |V> 3 ), IV'o), 
\ip'o))- -A- summary of this information is given in Fig. ([5]) 
of appendix A. 

In the following, we will classify the RG equation of 
the regions where each of this states represent the ground 
state. 
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A. Region (A): eo is the ground state. 

In this region the effective Hamiltonian in the first or- 
der correction is similar to the initial one, i.e, 



J' 



N 



-D'(a*aV +1 -aVa? +1 ) 



(3) 



where J', A' and D 1 are the renormalized coupling 
constants. The new renormalized coupling constants are 
found to be functions of the original ones given by the 
following equations, 
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»8 




J' = J(-) 2 (l + D 2 ), D' = D 

A f A + g 2 
1 + D 2[ 4 ' ' 



(4) 



A' 



FIG. 3: (color online) DC order (solid lines) and Staggered 
Magnetization (dotted lines) versus D. 



B. Region (B): 63 is the ground state. 

In this region the effective Hamiltonian to the first or- 
der corrections leads to the ferromagnetic Ising model 



H eff 



1 



N/3 



where 



A' = JA, J > 0, A < 0. 



III. PHASE DIAGRAM 

A. Region (A) 

For simplicity we have separated this region into posi- 
tive anisotropy and negative anisotropy sectors. 

• A > 

In the positive anisotropy sector the RG equations 
show that the J coupling, representing the energy 
scale, approaching zero by iterating RG procedure. 
Thus, at the zero temperature, the quantum phase 
transition is the result of competition between the 
anisotropy (A) and the DM coupling constant (D). 
In the region of planar anisotropy < A < 1, 
the symmetric interactions (D = 0) is known not 
to support any kind of long range order and the 
ground state is the so called spin-fluid (SF) state. 
Increasing the amount of anisotropy is necessary 
to stabilize the spin alignment. For A > 1 the 
ground state is the Neel ordered state. In the case 



of D ^ 0, the anisotropy constant (A) and antisym- 
metric (DM) coupling are in competition with each 
other. The latter thus destroys the ordering ten- 
dency of the former and defers creating of Neel or- 
der. Our RG equations show that the phase bound- 
ary between the SF and Neel phases which depends 
on the DM coupling is A c = y/l + D 2 (see Fig.©) 
which agrees with the phase boundary reported in 
ReflE This critical line coincide with boundary 
line which obtained by classical approximation (see 
appendix C). The RG equations (Eq.(H])) express 
the DM coupling dose not flow under RG transfor- 
mations, and the anisotropy coupling goes to zero 
(A — > 0) in SF phase while it scales to infinity 
(A — > 00) in the Neel phase. 

We have linearized the RG flow at the critical line 
A c = vl + D 2 and found one relevant and one 
marginal directions. The eigenvalues of the matrix 
of linearized flow are Ai = |, A2 = 1. The cor- 
responding eigenvectors in the |A, D) coordinates 
are |A X ) = |1,0), |A 2 ) = \^ m ,\). The marginal 
direction corresponds to the tangent line of the crit- 
ical line and the relevant direction shows the direc- 
tion of anisotropy 's flow (Fig.@). 

However we have found the boundary of the SF- 
Neel transition by calculating the staggered mag- 
netization Sm (sec appendix B) in the z-direction 
as an order parameter (Fig. ([3]) and Fig.(|4])), 



(-!)< 



(5) 



i=l 



Sm is zero in the SF phase and has a nonzero value 
in the Neel phase. Thus the staggered magnetiza- 
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FIG. 4: (color online) DC order (solid lines) and Staggered 
Magnetization (dotted lines) versus A. 



tion is the proper order parameter to represent the 
SF-Neel transition. We have plotted Sm versus D 
and versus A in Fig. ([3]) and Fig. ([4]), respectively. In 
Fig. ^ it is obvious that the staggered magnetiza- 
tion goes to zero continuously at the critical value 
of D c which shows the destruction of Neel order. 
The critical value D c where the staggered magne- 
tization vanishes above it, increases by increasing 
of anisotropy, this means SF-Neel transition point 
(A c ) depends on the DM interaction. It is seen 
in Fig. (|3]) that the staggered magnetization is zero 
for A < A c (SF phase) and has a nonzero value 
for A > A c (Neel phase), while enhancing of DM 
coupling defer the creation of Neel order. 

Moreover, to study the influence of DM coupling 
we have calculated the chiral orde r 26 i 27 (C^) in the 
z direction as the helical magnetization in one 
dimension 28 which is created by DM interaction. 

1 N 1 

i=i 

Unfortunately the chiral order is not a self similar 
operator under RG transformations. In fact an XX 
term of Hamiltonian (of of +1 + of of +1 ) shows up 
in the renormalized chiral order under RG. Thus, 
we have to calculate the sum of chiral and XX 
term which we have represented it by DC in the 
following equation, 



1 N 1 

i=l 



erf a 



t) 



erf o; 



i+ 



l)>- 



(6) 



In Fig.© and Fig.® the DC order has plotted 
versus D and A respectively, for different values of 
A and D. The figures manifest that the DC order 
enhances by increasing of DM coupling and reduces 
with increasing of anisotropy parameter. 

We have also calculated the critical exponents at 
the critical line (A c = \/l + D 2 ). In this respect, 
we have obtained the dynamical exponent, the ex- 
ponent of order parameter and the diverging ex- 
ponent of the correlation length. This corresponds 
to reach the critical point from the Neel phase by 
approaching A — > A c . The dynamical exponent is 
z ~ 0.73, the staggered magnetization goes to zero 
like S M ~ |A - A c f with (3 ~ 1.15. The correla- 
tion length diverges £ ~ |A — A c |~ 1 ' with exponent 
v ~ 2.15. The remarkable result of these exponents 
is the independence of their values on the D value 
and equality of them with the correspondening ones 
in the XX Z model. The detail of this calculation 
is similar to what is presented in Refi^. 

• A < 

In this sector, the effective Hamiltonian is simi- 
lar to the positive anisotropy case with the same 
coupling constants. For —yl + D 2 < A < the 
ground state is the spin-fluid phase and decreas- 
ing the anisotropy causes the ground state of the 
three site Hamiltonian changes by level crossing at 
A = — yl + D 2 where the RG equations should be 
reconstructed. However, the remarkable result is 
that the level crossing line which got by three site 
Hamiltonian coincides the critical line of this model 
in the thermodynamic system 2 ^. The RG equations 
express the DM coupling dose not flow under RG 
transformations, and the anisotropy coupling goes 
to zero (A — > 0). Thus the A = line is the po- 
sition of stable fixed points where the RG flow is 
freezed. 



B. Region (B) 

As we pointed out in sec[H]Tj the original Hamiltonian 
is mapped to the ferromagnetic Ising model. Ising model 
remains unchanged under RG as the stable fixed point 
and its properties are well known. We call this region as 
the ferromagnetic Ising phase. 



IV. SUMMERY AND CONCLUSIONS 

We have applied the RG transformation to obtain the 
phase diagram, staggered magnetization and helical mag- 
netization of XX Z model with DM interaction. In the 
positive anisotropy region, tuning the anisotropy cou- 
pling makes the system to fall into different phases, i.e 
Neel phase with nonzero order parameter and spin-fluid 
one with vanishing order parameter as characterized by 
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FIG. 5: (color online) The ground state eigenvalues as a func- 
tion of anisotropy and DM coupling. The thick long dashed 
line which shows the border line of region (A) and region (B) 
are given by A = — \/l + D 2 . 
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V. APPENDIX 



A. The block Hamiltonian of three sites, its 
eigenvectors and eigenvalues 



We have considered the three-site block (Fig. dT])) with 
the following Hamiltonian 



;, b J , j j 
h I = 4 { a lJ a 2J 



2J U 3J 



the staggered magnetization. The RG equations state 
that the system has fixed points at A = 0, A = oo and 
A c = \/T+ D 2 . The fixed points A = and A = oo 
are attractive and correspond to the Spin-Fluid and Neel 
phases, respectively. The fixed points at A c = vl + D 2 
are repulsive and correspond to the critical points of this 
model, in the other word it is the critical line of this 
Hamiltonian. However, in the negative anisotropy re- 
gion, the level crossing line A = —y/l + D 2 which is ob- 
tained by three site block Hamiltonian eigenvalues, is the 
critical line of infinite size system and separates the fer- 
romagnetic and spin-fluid phases. Unfortunately we can 
not calculate the chiral order explicitly by RG method. 
To survey the influence of DM interaction and helical 
magnetization, the numerical Lanczos computation is in 
progress. 
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The inter-block (H BB ) and intra-block (H B ) Hamil- 
tonian for the three sites decomposition are 

N/3 

i 

N/3 

i 

+A(of iJ af iJ + aljolj) 

+ d Ki4,i - + </4/ " °t,i°h) ■ 

where <r" 7 refers to the cv-component of the Pauli ma- 
trix at site j of the block labeled by /. The exact treat- 
ment of this Hamiltonian leads to four distinct eigen- 
values which are doubly degenerate. The ground, first, 
second and third excited state energies have the following 
expressions in terms of the coupling constants. 
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-A)(l 
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f £ 2 ) 


VMq - 


A)(W 


-L» 2 ) i 



2p 2 + 1)| ||t) - (1 - iD)(A + q)\ HI) - 2[2iD + (D 2 - 1)]| |||) 
2{D 2 + 1)| |tt) - (1 - iD)(A + q)\ HI) - 2[2iD + (D 2 - 1)]| TTi>) 



eo = --(A + g), 
2p 2 + 1)1 lit) - (1 - iD)(A - g)| HI) - 2[2iD + (D 2 
2(D 2 + 1)| |TT> - (1 - iD)(A - q)\ TIT) - 2[2iD + (D 2 - 



-1)]|TU> 
-1)]|TU)) 



eo 



I 02 



W2) = 



1 



V2(l + ^ 2 ) 
1 



V2(l + D 2 ) 



[2iD + (D 2 - \)]\ TII) + (^ 2 -l)| III) 

[2iD + (D 2 - \)]\ nj) + (D 2 -i)| ITT) 



e 2 =0, 

w = im) . w> = uu>, 

ea = £(A), 



(7) 



(8) 



(9) 



r 



where g is g = ^A 2 + 8(1 + L> 2 ). 

I T) and I I) are the eigenstates of a z . In Fig.© we 
have presented the different regions where the specified 
state is the ground state of the block Hamiltonian. In 
the region (A) the projection operator is 

Po = I Wo| + I Wol- 



The Pauli matrices in the effective Hilbert space have the 
following transformations 



P 1 n x P 1 



P 1 ^ P 1 
r a lJ r 



Do<\) 



P i n x pi _ 4(J) +1) IX j j 

--(Da, -a,) , P offw P =-__ ffj 

zq q 



P / fT W P 1 



-{Da j +a j) 



In the region (B) (A < —\J\ + D2) the projection op- 
erator is 



Po = \it)(H + I Wal- 

, u a TT . 1U L W = o - ^4^0=0 , py^ = o 

and the Pami matrices m the enective Hubert space have , T T , , , ., 

1 pi _2 pi pi _ pi _z pi _/* 

the following transformations ^0 ^u^o - *o °3,J-M> - ^0 cr 2,i^o — a i- 
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B. Order Parameter and Chiral Order 

1. Staggered magnetization 

Generally, any correlation function can be calculated 
in the QRG scheme. In this approach, the correlation 
function at each iteration of RG is connected to its value 
after an RG iteration. This will be continued to reach a 
controllable fixed point where we can obtain the value of 
the correlation function. The staggered magnetization in 
a direction can be written 

1 N f-lV 

s m = nY,(°\ CTf|0> ' (10) 

i 

where erf is the Pauli matrix in the ith site and \G) is the 
ground state of chain. The ground state of the renormal- 
ized chain is related to the ground state of the original 
one by the transformation, Pq\0') = \0). 

1 N f-lV 

i 

This leads to the staggered configuration in the renormal- 
ized chain. The staggered magnetization in z direction is 
obtained 

i N r-iv 

S% = ^£<0|^f-of|0) 

i=l 

= g Z E [{0'|P / (-<x+4/-4/) P o|0') 

3 1=1 

- (o'|p /+1 + - </+i)n /+1 |oO 

2A + g i v\ n ,,(-i)* ^m'\ A* nn 

= o — iv z_^( I o i\ i = — t^mA 11 ) 

6q T i=i 1 6 

(n) 

where S M is the staggered magnetization at the nth step 
of QRG and is defined by 7 = (2 A + q)/q. 

This process will be iterated many times by replac- 
ing with 7W. The expression for 7'™' is similar to 
<y(°) where the coupling constants should be replaced by 
the renormalized ones at the corresponding RG iteration 
(n) . The result of this calculation has been presented in 
Fig.© and Fig.®. 

2. Chiral Order 

The chiral order which is the proper function to detect 
the helical magnetization in the systems can be written 



As we mentioned in the section III, the XX term of 
the Hamiltonian shows up to the chiral order under RG. 
The XX term order is written 



Dxx = jj £ 4<«<+i + tM+i)>- (13) 
i=i 

In this case the calculating of the chiral order being 
elaborate, because of the unknown effect of the XX term 
on the ground state of system at fixed point (A = 00). To 
simplification of the calculation, we transform the XX Z 
with DM interaction Hamiltonian (Eq.®) to the Ising 
model with DM interaction 28 (IDM) by implement a non- 
local transformation which shows a rap rotation about 
the z axis at site n where ip = arctan(— -^). We have 
calculated the chiral order fEq.([T2"])) of IDM in Rcf.[ 28 ]. 
By implement the inverse of transformation, the chiral 
order in IDM model transforms to the DC order where 
introduced in Eq.©. The DC order has been plotted in 
Fig.© and Fig.® versus D and A. 



C. Classical Approximation 

In the classical approximation the spins are considered 
as classical vectors which form the spiral structure with 
a pitch angle if between neighboring spins and canted 
angle 9 

a* = cos(n(p) sin , a\ = sm(rup) sin , o z n = cos 8, 

The classical energy per site for the XXZ with DM 
interaction Hamiltonian (Eq.|T])) is 

E J 

-£ = -7 [( cos V + D sin <f) sin2 + A cos2 O] ■ 

The minimization of classical energy with respect to 
the angles tp and 8 shows that there are two different 
regions. (I) A > vl + D 2 , the minimum of energy is 
obtained by arbitrary 8 and ip — arctan(D) which show 
the spins projection on z axis is nonzero and spins have 
the helical structure (see Fig|5|) in the xy plain. In this 
region the minimum classical energy is 

F 1 J 

= -(v / l + £ )2 sin 2 6' + Acos 2 8). (14) 

(II) A < \J\ + D 2 , the energy is minimized by A = 
cos ip + D sin tp and arbitrary 9 or arbitrary <p and 9 = ? , 
which correspond respectively to the configurations with 
nonzero value of spins projection on z-axis with helical 
structure of spins projection in the xy-plain and disor- 
der configuration. In this region the minimum classical 
energy is 
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FIG. 6: (color online) A classical picture of spin orientation in 
the xy plain where the angle between neibouring spins depend 
on the D value. 



N 
3=1 



G 3 



J-l 

ctj — mt&TT 1 {D), 

m— 1 



H = e- lU He iU , 



(16) 



El? J , „ . . J . 

— 9L. = — (costp + Dsmip) = — A. 



(15) 



One can see from Eq. (|T4| and Eq. (|15|) that the tran- 
sition between phase (I) and (II) takes place at A = 

vtttp~. 



VI. CANONICAL TRANSFORMATION 



which gives 



H = 



j yi + tp 



VT+TP' 



(17) 

The f/(l) x Z2 symmetry of initial Hamiltonian sur- 
vive in the transformed Hamiltonian too, but at A c = 
±Vl + D 2 the L^(l) x Z2 symmetry breaks to the local 
SU(2) symmetry. 



The Hamiltonian of XX Z model with DM interaction 
(Eq.flTJ) has the global U(l) x Z 2 symmetry. This Hamil- 
tonian is mapped to the well known XX Z chain via a 
canonical transformation 25-29 , 
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